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I. INTRODUCTION 



Leading three-loop logarithm cubed and logarithm squared radiative-recoil contributions 
to hyperfine splitting (HFS) in muonium were calculated long time ago (see, e.g., reviews in 
[1112]). Recently we started calculation of all single-logarithmic and nonlogarithmic radiative- 
recoil correction (see review in |3| ) . Below we consider single-logarithmic and nonlogarithmic 
radiative-recoil corrections due to insertions of electron and muon polarization operators in 
the radiative photon in Figs. [T} [2j 




FIG. 1: Electron polarization insertions 
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FIG. 2: Muon polarization insertions 

Three-loop diagrams in Figs. [l| [2] can be obtained from the diagrams with two-photon 
exchanges by insertion of radiative corrections in Fig. |3| The two-photon diagrams in Fig. [3] 
produce the leading radiative-recoil correction when the loop momentum is much larger than 
the electron mass, and insertion of radiative correction can only increase the integration 
momentum. Therefore, calculating the diagrams in Figs. [TJ|2]we may forget about external 
virtualities and calculate matrix elements in the scattering regime between the free electron 
and muon spinors. To turn the matrix element into contribution to HFS we multiply the 
scattering matrix element by the Coulomb- Schrodinger bound state wave function at the 
origin squared, and calculate difference between spin one and spin zero states. We use 
the Feynman gauge to obtain matrix elements of the gauge invariant sets of diagrams in 
Figs. [T| [2} Each of the dia grams in Figs. [TJ |2] contains polarization operator insertion in 
one of the radiative photon lines. We account for this insertion using the massive photon 
propagator for radiative photons (but not for exchanged photons) with the photon mass 
squared = 4m^/(l — v'^) or = 4M^/(1 — v"^), where m and M are the electron and 
muon masses, respectively. Insertion of the polarization operator in the radiative photon 
line is accompanied by an additional integration over velocity v with the weight 
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FIG. 3: Two-photon exchanges 

Besides single-logarithmic and nonlogarithmic contributions the diagrams Fig. [T] generate 
also well known much larger nonrecoil and logarithm squared recoil contributions [1]. We 
will calculate the contributions of the diagrams Figs. [T| |2] with linear accuracy in the small 
electron-muon mass ratio m/M. In particular we will reproduce nonrecoil and logarithm 
squared recoil contributions what will serve as an additional check of our new results. The 
paper is organized as follows. In Section II we describe calculations of the diagrams with 
the electron polarizations in Fig. [T| and Section III deals with the diagrams with the muon 
polarizations in Fig. [2j All results are collected in the last section. 




II. ELECTRON POLARIZATION OPERATOR 

A. Calculation of the Mass Operator Contribution 

We calculate matrix elements of each of the diagrams in Fig. [T] separately. Respective 
integrals can be obtained by inserting radiative corrections in the expression for the con- 
tribution of the skeleton diagrams in Fig. [3] to HFS. Contribution of the diagram with the 
self-energy insertion in Fig. [ijhas the form (compare P]) 
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where dimensionless energy Aes is connected with the energy shift AEy. by the relation- 
ship AEj] = {a'^{Za)/n^){m/M)EFAeT,^, fi = m/{2M), k is the dimensionless Minkowski 
exchanged momentum, Aesi and Aes2 are the integrals corresponding to two terms in the 
square brackets, and 
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The integration over v in Eq. ^ take s care the electron loop that is included in the integrand 
via the finite mass A = a/4/(1 — v'^) of the radiative photon. We are going to calculate two 
leading terms in the expansion of the dimensionless energy Ae^ in Eq. ^ over the small 
parameter /i (the term of order l/yu and the term independent of /i). 



^ The Fermi energy is defined as Ep — {8/3){Za)'^{m/M)m. 
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1. Nonrelativistic Contribution 



As a first step of calculations we obtain tlie leading term of order from the expres- 
sion for Aes in Eq. This term gives the leading nonrecoil contribution to HPS (recall 
AE-£ = {a'^{Za)/TT^)2fiAeY,) and arises, as all nonrecoil contributions, in the external field 
approximation. To extract the nonrecoil contribution from the expression in Eq. ^ we take 
the residue at the muon pole what with linear accuracy in fi reduces to substitution 



2k^ 



- fi-^k^ 

Then we obtain the nonrelativistic for the muon contribution in the form 
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This last integral contains both recoil and nonrecoil contributions. The recoil contribution 
will be treated on equal grounds with other recoil contributions considered below. Integrat- 
ing in Eq. ([s]) over momentum and expanding the result over the small parameter fj, we 
obtain 
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2. fi - and c - Integrals 

We return now to calculation of the first two terms in the expansion of the contributions 
to HPS in Eq. Q over fx. An attempt to calculate the integral in Eq. ^ with the help of 
the Peynman parameters leads to the integrands that do not admit expansion over small 
parameter /i before integration. Therefore we use another approach to calculation of the 
integral in Eq. ^ (as well as to calculation of other integrals of this type below), and directly 
integrate over the exchanged momentum k in four-dimensional polar coordinates. After the 
Wick rotation and integration over angles (and omitting some higher order terms in /i) we 
obtain integral representation for the contributions Ae^i and Aes2 defined in Eq. ([2]). The 
integral representation for Ae^i (compare |5j) has the form 
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and the integral representation for Aes2 (compare has the form 
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Each integrand in Eq. ([t]) and Eq. ([s]) is a sum of /i-dependent and /i-independent integrals, 
and we would like to calculate them separately. While both integrals are convergent, naive 
separation of /x-dependent and /^-independent terms in the integrands sometimes leads to 
integrals that are ultravioletly divergent at large integration momenta k. For example, 
integral over k of /i-dependent term in the integrand in Eq. ([T]) converges at large integration 
momenta, while respective integral in Eq. ([s]) diverges. This divergence arises because fi- 
dependent terms in the integrand go to /^-independent constants at high momenta. In 
such cases we redefine the /i-dependent terms in the integrand by subtracting the leading 
asymptotic constant, and add this constant to the /i-independent terms in the integrand. A 
universal recipe for such restructuring of the integrands in Eq. ([T]), Eq. ([s]) is described by 
the substitutions 
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In the first case we subtract 1, in the second case we subtract 1/2. In both cases we add the 
terms corresponding to these constants to the /i-independent terms in the integrands. After 
this restructuring (when needed) we write the integrals Eq. ([T]) and Eq. ([s]) as sums of what 
we will call /i- and c-integrals 



(10) 



Let us consider calculation of /t-integrals. The integral over k in Eq. Q converges, the 
integrand does not require any restructuring, and the /i-integral has the form 
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This integral is of order /i, and therefore with our accuracy it does not give any contribution 
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Next we calculate the /i-integral arising from the integral in Eq. ([8 
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Unlike the case of Ae^i this time the naive integral over k of the /x-dependent terms in 
Eq. ^ diverges at large k, and we restructured the integrand according to Eq. Besides 
the recoil contribution the integral in Eq. ( [I3| ) contains also nonrelativistic contribution 
Ae^(A^-R) that we calculated separately in Eq. This nonrelativistic contribution is 

generated by the leading l/{fik) term in the expansion over small fik of the expression in 



the square brackets in Eq. ( 13 ). It coincides with the contribution generated by the first term 
in the square brackets in Eq. To avoid double counting we subtract this contribution 
from the integrand in Eq. (13) by the substitution in the integrand 



(14) 



This substitution gives a universal recipe for subtraction of the nonrecoil corrections in all 
/i-integrals to be considered below, and it is necessary to emphasize that it is needed only in 
the first of two typical structures with square roots in Eq. (|9| that arise in the expressions 
for /^-integrals. The leading term in the small fik expansion of the second structure is 
nonsingular and does not generate nonrecoil contribution. 

Finally, the second /^-integral for the mass operator insertion in the electron line has the 
form 
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Other integrals of this type will arise below in calculations of other contributions to HPS. 
To extract the first term in the expansion of integrals of this type over the small parameter 
fi we introduce an auxiliary parameter a, such that 1 -C a -C l//i (see, e.g., [HE]). Then we 
separate the large and small integration momenta regions with the help of this parameter 
a and use different approximations in different regions. In the region of small integration 
momenta < < cr we expand the integrand over fik <^ 1 and obtain 



AeS<^ - A£;(r^(iV^) -^lnV + glnV+ (^-n' + ^) lna + 1.3220. (16) 

In the region of large integration momenta k < a we expand the integrand over the small 
1/fc -C 1 and obtain 
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In the intermediate region k a both approximations fik <^ 1 and l//c -C 1 are valid 
simultaneously, and all dependence on the auxiliary parameter a cancels in the sum of the 



contributions in Eq. (16) and Eq. (17) 
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Let us turn to calculation of the c-integrals as defined in Eq. (10), Eq. ([t]), and Eq. 
We easily perform momentum integration in the integral for Ae^{ 
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As in the case of the /i-integrals above we want to avoid double counting, and subtract the 
respective nonrelativistic contribution already accounted for in Eq. ([6]). This nonrelativistic 
contribution has the form 



pl fX pi 

Ae^^liNR) = 3 dx dy - 
Jo Jo Jo 1 



dv 



v'[ 1 



v^\ hi{x,y) 
3 / ai 



(20) 



Now we see that due to the identity tan ^(1/ai) = 7r/2 — tan ^ ai subtraction of nonrela- 
tivistic contribution from the integral in Eq. (|T9|) reduces to substitution 
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This is a universal rule for subtraction of nonrelativistic contributions in all c-integrals 
considered below (compare [S]). Finally, the first c-integral is 
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Next we turn to the second c-integral defined in Eq. (10) and Eq. ([8j), and start with the 
momentum integration (compare [3]) 



Ae 



S2 



dx / dy 



dv 



1 



V 

~3 



h2ix,y) 



dk^ 



(23) 



X 



dx / dy 



dv 



h2ix,y) 



2 1 1 3 1 + a? 1 a2 1 + a? 
tan" — + -ln ^ + -- — \n ^ 



ai 4 



a 



4 4 



a 



7 



Subtraction of the nonrecoil contribution is needed to avoid double counting, and it is done 
with the help of the universal rule in Eq. (21). Then we obtain 
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Next we collect all contributions in Eq. Q, Eq. Eq. (|T8|, Eq. and Eq. (|24|), 

and obtain the total result for the contribution of the diagram with the self-energy insertion 
in Fig. [T] 
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B. Calculation of the Spanning Photon Contribution 

Contribution of the spanning photon diagrams in Fig. [T] is obtained from the contribution 
of the skeleton diagrams in Fig. [3] by insertion of the radiative photon with the polarization 
bubble, and is described by the expression (compare |1]) 
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We simplify this expression using the identities 
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TABLE I: Spanning Photon Contributions 
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where the contributions Ae=i and Ae=2 correspond to the expressions in the first and second 



square brackets on the right hand side in Eq. (29) 



For calculation of the of the integrals in Eq. (29) we use the same tricks as in the case of 
the mass operator contribution. We skip the details of calculations and collect intermediate 
results in Table 1. Finally we obtain the contribution of the diagram with the spanning 
photon insertion in Fig. [T] in the form 

A6H = ^ + \ ln^(2/.) - ^ ln^(2;.) + + |) ln(2/i) - 0.6314. (30) 



C. Calculation of the Vertex Contribution 

Contribution of the diagram with the vertex insertion in Fig. [T] is obtained from the 
contribution of the skeleton diagrams in Fig. [3] by substitution of the vertex function instead 
of one of the skeleton vertices. We have derived a convenient expression for the one-loop 
vertex function with a massive photon 
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TABLE II: Vertex Contributions 
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where Aq{x) = + A^(1 — x), and the terms A^*-* correspond to the five terms in the braces 
in Eq. 

This is essentially the same expression as the one in But unlike the respective ex- 
pression in [1], where the photon mass merely served as a regularization parameter and was 
preserved only when necessary, here we have restored full dependence on the finite photon 
mass A that effectively describes polarization operator insertion. One can show that the 
gauge invariant anomalous magnetic moment does not generate radiative-recoil corrections 
(see, e.g., Ej). Therefore, the anomalous magnetic moment and some other terms that 
do not contribute to HFS are omitted in Eq. (31). 

We plugged in the vertex in Eq. (31) in the skeleton expression for the contribution to 
HFS, and obtained five integrals corresponding to the five terms on the right hand side in 
Eq. (31). Calculation of these integrals goes along the same lines as in the case of the mass 
operator, discussed in detail above. We collected all intermediate results in Table II. Total 
contribution of the diagram with the vertex insertion in Fig. [l] turns out to be 
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Now we collect all contributions in Eq. (25), Eq. ipOh, and Eq. (32) generated by the 



diagrams with electron polarization insertions in Fig. Ol and obtain 



1 5335 3 / 53\ 

Ae^'^) = Aes + 2AeA + Ae= = + - ln'(2/i) + -tt' + — ln(2^) + 7.0807. (33) 

4 V 6 7 

The first term on the right hand side is the well known nonrecoil contribution to HFS 
[8], the second term is the leading logarithm squared contribution obtained in [9], and the 
single-logarithmic and constant terms are the subject of this work. 



III. MUON POLARIZATION OPERATOR 

Insertion of the muon polarization operator in the radiative photons in the diagrams in 
Fig. |2] lifts characteristic integration momenta to the scale of the muon mass. Hence, these 
diagrams do not generate nonrecoil contributions to HFS, all of which originate from the 
region of nonrelativistic muon momenta. Moreover, due to high characteristic momenta 
these diagrams even do not generate logarithmic in the mass ratio recoil contributions that 
originate from the wide integration region between the electron and muon mass. As a result 
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leading recoil contributions of the diagrams in Fig. [2] are pure numbers, and their calculation 
is significantly simpler than in the case of the electron polarization insertions in Fig. [TJ 

Like in the previous section, insertion of the muon polarization operator in the diagrams 
in Fig. [2] is accounted for by introduction of a photon mass, followed by an additional 
integration over the velocity with the weight from Eq. ([T]). For muon polarization the 
effective photon mass in the integrals is large, in dimensionful units = 4M^/ (1— f ^), and it 
determines characteristic momenta in all integrals. We can obtain the expressions for energy 
shift due to muon polarization from the formulae for the respective electron polarization 
contributions above by rescaling the dimensionless integration momenta k ^ k/fi (recall that 
fi = m/{2M)). In addition we should adjust the expression for the photon mass, in terms of 
the rescaled integration momenta measured in units of 2M it is = l/(/i^(l — f^)). After 
these substitutions the expressions for the dimensionless contributions to HFS that are due 
to electron polarization go into expressions for the contributions due to muon polarization. 

a. Mass operator contribution. To obtain an explicit expression for the diagrams with 
the mass operator insertions in Fig. |2] we rescale the integration momentum k ^ k/jj, m. 
Eq. ([7| and Eq. ([s]), and redefine the photon mass squared as A^ = l/(/i^(l — f^)). The 
auxiliary functions used in Eq. ([T]) and Eq. ([s]) are defined in Eq. (|3]). After rescaling they 
simplify 
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All dependence on fi becomes explicit after these manipulations. It turns out that Ae^i van- 
ishes together with /i. The total leading recoil contribution generated by the diagrams with 
mass operator insertions in Fig. [2] coincides with Ae-S2. Its calculation is straightforward, 
and we obtain 
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h. Spanning photon contribution. We obtain an expression for the spanning photon 
contribution with the muon polarization insertion in Fig. |2] by rescaling the integration 



function A(x,y) in Eq. (27) simplifies 



momentum and the photon mass in Eq. (|29|). Under these transformations the auxiliary 
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After rescaling only the first and third terms in the first square bracket on the right hand 
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side in Eq. (29) produce nonvanishing with contributions, and we obtain 
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^(^Vl + k^-k^ - ^(^kVlT¥ - P - =0.3105. 
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c. Vertex contribution. Rescaling the integration momentum and the photon mass in 
the expressions corresponding to the five terms A^*'' we obtain an exphcit expression for the 
vertex diagram contribution in Fig. [2] to HFS. It is easy to see that after rescahng all terms 
in the expression for the vertex function in Eq. (31) are suppressed by at least one power 
of fi in comparison with the first term. This means that only the first term generates the 
leading recoil correction in the case of muon polarization insertion. Explicitly, the leading 
recoil correction generated by the vertex insertions in Fig. |2] has the form 



dx / dy 



dvv^ 



V 

~3 



' {x-y){l-2y)+y{l-y) 

kM^-y){l-v^) + l-x ^ ' 

1 



-0. 8738. 



Collecting all leading radiative-recoil corrections in Eq. (35), Eq. (37), and Eq. (38) 
corresponding to the diagrams with muon polarization insertions in Fig. [2| we obtain 



Ae^^) = Ace + 2AeA + Aes = -1.3042. 



(39) 



IV. CONCLUSIONS 



Restoring the overall dimensional factor in Eq. (33) and throwing away nonrecoil and 



logarithm squared terms known earlier, we obtain single-logarithmic and nonlogarithmic 
contributions to HFS generated by the diagrams with one-loop electron polarization inser- 
tions in Fig. [T] 
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The radiative-recoil contribution generated by the diagrams with one-loop muon polarization 
insertions in Fig. |2]is nonlogarithmic. We obtain it restoring the overall dimensional factor 



in Eq. (39) 



-1.3042 ^ !^ -—Ep. 
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(41) 



The total contribution to HFS obtained above can be written as (Z = 1 in muonium) 



AE = AE^^) + AE^^"^ 



TT 
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m 
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(42) 



Currently the theoretical accuracy of HFS in muonium is about 70 Hz. A realistic goal is 
to reduce this uncertainty below 10 Hz (see a more detailed discussion in [H 12]). The result 
in Eq. (42) together with other three-loop radiative- recoil results in [TUl [TTl |T2] makes this 
goal closer. 
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